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Abstract 

In this paper a new version of the DHN (Dashcn-Hasslacher-Neveu) formula, which is used to 
compute the one-foop order kink mass correction in (l+l)-dimensional scalar field theory models, is 
constructed. The new expression is written in terms of the spectral data coming from the supersym- 
metric partner operator of the second-order small kink fluctuation operator and allows us to compute 
the kink mass quantum shift in new models for which this calculation was out of reach by means of 
the old formula. 

PACS: 11.15.Kc; 11.27.+d; ll.lO.Gh; ll.30.Pb 

1 Introduction 

In this work we shall bring together two research subjects in theoretical physics grown respectively 
in the late seventies/early eighties and developed along the past forty years. The first topic is the 
study of quantum fluctuations of kinks and other topological defects. Dashen, Hasslacher, and Neveu 
(DHN) [ ] started the investigation of this problem within the framework of the ^-expansion in Quantum 
Field Theory with the motivation of understanding hadrons as quantum descendants of solitary (or even 
solitonic) non-linear waves. Very soon after, several authors shaped this field of research almost definite 
both in Russia, Faddeev and Korepin [2], and the United States, Coleman, Jackiw, Goldstone, Christ 
and Lee [3, 4, 5, 6]. The second theme is supersymmetric quantum mechanics. Introduced by Witten 
[7] in order to elucidate a mechanism of spontaneous supersymmetry breaking, the idea triggered an 
infiationary process of exploring physical supersymmetric quantum mechanical systems, see the review 
of Cooper et al. [o] and the books of Freund [y] and Junker [iU] to find full information about the basic 
ideas in SUSY quantum mechanics. 

In references [11, 12] DHN performed a very detailed analysis of the kink mass shift issue in the (1 
+ l)-dimensional sine-Gordon and (p'^ scalar field theory models. In particular, the authors succeeded 
in computing the one-loop order mass shifts for the soliton and kink non-linear waves existing in these 
models. They recognized that the kink mass shift induced by the quantum fluctuations must be due 
to three contributions: 1) the kink zero-point energy collecting the energy of the kink ground state, 
where all the fluctuation modes are unoccupied, 2) the analogous vacuum zero-point energy that must 
be subtracted, and 3) the energy induced by the one- loop renormalization mass counter-term on the 
kink background (measured with respect to the same effect on the vacuum). Note that the first two 
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contributions require the computation of the difference between the traces of the kink and vacuum 
fluctuation operators, which are respectively Schrodinger and Helmholtz differential operators. The final 
expression encompassing all these contributions is known as the DHN (Dashen-Hasslacher-Neveu) formula 
and expresses the kink mass quantum shift in terms of the spectral data -bound state eigenvalues and 
scattering wave phase shifts- of the second-order small kink fluctuation operator. A remarkable fact to 
be stressed is that the kink fluctuation operators governing the 0^ and sine-Gordon kink fluctuations 
are reflectionless Poschl- Teller type Schrodinger operators. In these cases, both the kink and vacuum 
fluctuation operators exhibit a half-bound state, a singlet state with eigenvalue just at the threshold of 
the continuous spectrum. A general study of the (1+1) dimensional scalar field theory models enjoying 
this type of kink fluctuation operators can be found in the references [6, 1-3, 14, 15]. 

As early as in 1978 Olive and Witten [16] unveiled the character of BPS states of the kinks (and 
other topological defects) in field theories with extended supersymmetry and the connection with the 
central charges of SUSY algebra. This property makes the BPS kink mass very robust against quantum 
fluctuations. Nevertheless, by the end of the decade of the nineties the topic of quantum fluctuations of 
kinks underwent a comeback in the supersymmetric fleld theoretical framework. A thorough analysis of 
the DHN formula is given in [17, 18], where the authors uncover the mode number cut-off regularization 
procedure underlying this formula. In this regularization method the difference between the traces 
of the kink and vacuum fluctuation operator is determined by computing the difference between the 
contributions of the same number of eigenmodes in the kink and vacuum fluctuation operators. The main 
goal of these works, however, was the quantization of supersymmetric kinks, where the contributions of 
fluctuations of the bosonic and the fermionic flelds cancel each other, relieving the need of identifying 
the details of the spectrum of the kink and vacuum Hessian operators. Subtleties, however, come into 
play which are related to: the choice of different boundary conditions, e.g., periodic boundary conditions, 
Dirichlet, and/or Robin boundary conditions, etcetera, ii) diverse regularization methods, for instance, 
energy cutoff, mode number cutoff, higher-order derivatives, iii) the development of phase shift analysis, 
see references [19, 20, 21, 22, 23]. In the (/>^ and sine-Gordon models, the contributions of the half-bound 
states to the one-loop mass shift due to kink and vacuum fluctuation operators annihilate each other 
in the mode number cut-off regularization context. Thus, the original DHN formula, which is perfectly 
adapted to reflectionless kink fluctuation operators, works properly. This old formula, however, needs to 
be slightly modifled to deal with models whose kink fluctuation operator has a non-vanishing reflection 
coefficient. In this case, the spectrum of kink fluctuations does not embrace a half-bound state, unlike 
the vacuum fluctuation operator. The necessary modiflcation was unveiled from the one-dimensional 
Levinson theorem and described in reference [2 -J]. The outcome is a generalized DHN formula, which is 
the starting point in this work. 

The DHN formula, old or new, demands full information about the second-order small kink fluctuation 
operator spectrum: not only the eigenvalues but analytical knowledge of the eigenfunctions, which encode 
the phase shifts, is necessary because the spectral densities enter the DHN formula. This is an important 
handicap in the effectiveness of these formulas since the number of exactly solvable spectral problems is 
reduced. Although this paper is devoted to the exact computation of the kink mass quantum correction 
by employing the DHN formula, we do mention the existence of alternative methods that approximately 
compute the one-loop kink mass shift from the spectral zeta function regularization method, see [25, 26, 
27, 28]. The idea is to rely on the determination of the spectral zeta function of the kink fluctuation 
operator via the Mellin transform of the heat trace of the same operator and then use the asymptotic 
expansion of this latter spectral function to obtain information about the kink fluctuation asymptotics. 
During the last decade two of us and our colleagues benefited from the theoretical machinery available 
by applying these tools at our disposal to write the one-loop kink mass shift as a truncated series in 
the Seeley coefficients of the kink-Hessian- heat function in several (1 + l)-dimensional scalar field theory 
models [29, 30, 31] and Abelian gauge Higgs models [32, 33, 34]. 
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In this paper, however, we shall profit from the spectral information about the kink fluctuation 
operator obtained by using the techniques available in supersymmetric quantum mechanics [S]. In these 
corresponding systems, the spectrum (eigenvalues and eigenfunctions) of a given Schrodinger operator is 
related to the spectrum of a SUSY partner operator: the two operators are iso-spectral (except possible 
zero modes) and the eigenfunctions are intertwined by the action of the supercharges. The spectrum of 
stable kink fluctuations is non-negative. Since it is the spatial derivative of the kink solution, the ground 
state (a zero mode) is always known^. Knowledge of the ground state wave function allows us to factorize 
the Schrodinger operator in two first-order differential operators (this being the building blocks of the 
supercharges) where the "superpotential "is the logarithm of the zero mode. The SUSY partner second- 
order differential operator is obtained by multiplying the two first-order operators in reverse order. By 
construction the SUSY partner operators are, up to zero modes, "quasi "-isospectral but frequently the 
eigenfunctions of one of the operators are simpler than the eigenfunctions of the other. 

Our strategy in this work is to rewrite the DHN formula in terms of the spectral data of the new 
operator in the case where the eigenfunctions of the original operators are too complicated to use the 
associated spectral density in the analytical integrals. In the privileged cases of sine-Gordon, cj)'^, and 
the so-called parent models, see [15], the SUSY partner operators describing the kink fluctuations enjoy 
exactly solvable spectral problems. Moreover, all these Schrodinger operators belong to the hierarchy of 
reflectionless Poschl- Teller potentials. Starting from a given kink Hessian one can descend through all 
the members of the hierarchy to end in the free Helmholtz operator. The necessary spectral data in the 
DHN formula applied to these kinks can be obtained by this process from the vacuum spectral data, such 
that the computation of the kink mass quantum correction becomes almost trivial. 

There are other models with much more complicate kink solitary waves where this SUSY quantum 
mechanical structure simplifies (enables) the computation of kink mass shifts due to one-loop fluctuations 
out of reach from a direct approach. Usually, the recursive use of the above scheme allows us to write 
the kink mass shift in terms of the spectral data for a SUSY partner operator that lacks bound states. 
Finally, we shall construct families of field theory models whose kink fluctuation operators share a unique 
SUSY partner operator. Therefore, the mass quantum corrections for all the kinks in these models are 
identical, forming a kind of SUSY universality class. 

The organization of this paper is as follows: in section 2 we introduce the notation and general 
background on the subject: kink mass quantum corrections in (l+l)-dimensional scalar field theory 
models. In this section we also make explicit the relations between the spectral data of the second-order 
small kink fluctuation operator and its SUSY-partner operator. In subsection 2.2 we re-express the 
generalized DHN formula in terms of the spectral information of this new operator, while in subsection 
2.3 we accomplish the same task for any operator in the hierarchy of SUSY-partner operators. We apply 
this reformulation of the DHN formula to the sine-Gordon, and parent models, profusely studied in 
the literature. The novel approach unveils the exact results in a more direct way than the standard 
method. Section 3 is devoted to the application of this scheme to the computation of the kink mass 
quantum correction in more sophisticated scalar field theory models. The first one is characterized by 
a potential energy density U{(j)) = |(/)^ cos^ log |(/)|, whose topological kinks have been studied in [35]. 
We shall compute the kink mass quantum correction in this model by following the strategy explained 
above. The spectral problem of the SUSY partner operator offers much more analytical information about 
the eigenfunctions than the original one. In subsection 3.2 we introduce a model reconstructed from a 
zero mode following the procedure described in [6]. In this case, the intrincate second-order small kink 
fluctuation operator has (miraculously) a Poschl- Teller type SUSY-partner operator, which again allows 
us to obtain the kink mass quantum shift. Finally, in Section 4 we construct families of models whose 
second-order small kink fluctuation operators share the same SUSY-partner operator. We shall describe 
in sub-Section 4.1 a family of one-gap Schrodinger operators having a unique common supersymmetric 

'^This is due to the fact that the spatial translational symmetry of the system is spontaneously broken by the kink. 
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partner. In sub-Sections 4.2 and 4.3 a similar analysis will be performed on two-parametric families of 
two-gap Schrodinger operators exhibiting an interesting hierarchical structure. 



2 SUSY partner operators and kink mass quantum corrections 

In this section we shall introduce the theoretical background that will be used to deduce a new version of 
the Dashen-Hasslacher-Neveu (DHN) formula depending on the spectral information of the SUSY partner 
to the second-order kink fluctuation operator. In some models the new formula offers us the possibility of 
successfully computing the one-loop kink mass quantum correction even though it cannot be calculated 
directly from the standard formula. 



2.1 Field equations and topological kinks 



The action functional governing the dynamics in our (l+l)-dimensional relativistic one-scalar field the- 
oretical models is of the form: 



1 dtp 

2 9y; 



dtp 
dtf 



Here, tp = t and = y are local coordinates in M^'^, which is equipped with a metric tensor r/^jy = 
diag(l, —1) and tp{y^^) : M}'^ — M is a real scalar field. 

We shall work in a system of units where the speed of light is set to one, c = 1, but we shall keep the 
Planck constant h explicit because we shall search for one-loop quantum corrections, proportional to h, to 
the classical kink masses. In this system of units, the physical dimensions are: [h] = [S] = ML, [y^J = L, 
[tp] = M2L2, [U] = ML-^. The models that we shall consider are distinguished by different choices of 
the part of the potential energy density that is independent of the field spatial derivatives: U[tp{y^^y\. 
In all of them there will be two special parameters, rrid and 7^, to be fixed in each model, carrying the 
physical dimensions: [m^] = L^^ and [7^] = M^^L^^. We define the non-dimensional coordinates, fields 

and potential in terms of these parameters: Xn = mdyn, xq = t, xi = x, (p = ^d'^P, U{(p) = -^U{tp). The 
action S[tp] is also proportional to a non-dimensional action, that is, S[tp] = ^S[(p] with 



S[<P] 



dx' 



d^cpdf'cp-uicp) 



(1) 



where (p{x^) : M}''^ — )• R is a non-dimensional real scalar. The second-order non linear field equations are: 

6U 



d'^cp d'^cp _ 
dt"^ dx'^ 5(j) 

The non-dimensional potential energy functional E[tp] = ^E[4>\ reads^: 



(2) 



E\ 



dx e{x) 



dx 



1 

2\dx 



+ U{(P) 



(3) 



To guarantee the existence of topological kink solutions we assume that U{(p) is a non-negative twice- 
differentiable function, i.e., U{(p) £ C^(M) and U{(p) > 0, with a discrete set of zeroes Ai: 



M = {<P^^^ eR:Ui(P'^') = 0, j = l,...,iV} 



(4) 



We shall work with non-dimensional quantities, fields, and parameters in the sequel. To recover the correct fully 
dimensional magnitudes we shall use the relations written above in the text. 
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These time-independent and homogeneous configurations (4) are the simplest solutions of the field 



equations (2) and the absolute minima of the energy: E[(tU^] = 0. The smah vacuum fluctuations 



SfQ\t,x) = (j){t,x) — (py^ SiVe also solutions of the second-order equations (2) if the linearized second- 
order equations are satisfied: 



The general solution of the linear PDE (5) is the Fourier integral transform 

^fo\t^^) = ^J [c,(fc)e-^*-^^- + c*(fc)e— (6) 

such that the dispersion relation u^j{k) = k"^ + v'j holds. Canonical quantization transmutes the Fourier 

coefficients to annihilation and creation operators, and the Fock space built from the vacuum (j)y = 
(Oj|(/)(t, x)|Oj) is the space of the fundamental quanta of the system. 

The configuration space C = {0(to,x) G Maps (M\M) /£;[(/>] < oo} is the union of infinite disconnected 
pieces: card(A^) > 1. If this is the case, spatially extended field configuration solutions of finite energy 
can exist. In their center of mass they solve the second-order ODE 



(7) 



together with the asymptotic conditions 



lim (t){x) = G X , lim = . (8) 



Because is a discrete set the evolution in time from a given set of asymptotic conditions (8) to another 
set would cost infinite energy. In particular, when = j_ it 1 there are kink-shaped topological solutions 
to the ODE (7), whereas if j_ = j+, besides the constant solutions there might exist non-topological 
(bell-shaped ) kinks. The topological kinks verify the first-order ODE: 

f^=±^wW) . (9) 

One easily proves that the solutions of this Bogomolny-Prasad-Sommerfield equation are also solutions 
of the ODE (7). From (9) one infers that the BPS kinks are monotonically increasing functions asymp- 
totically connecting two consecutive minima, whereas the BPS anti-kinks are monotonically decreasing 
functions connecting the same pair of minima. 

A Lorentz transformation sends the static solution 4>k{x) to (j)K{t,x) = 4>K{-^f^^)- Thus, kinks 
are traveling waves that spontaneously breaks Lorentz invariance. Also, the spatial translations and 
refiections, x — >• x + xq, x i— )■ — x, leave the action invariant because there is no explicit dependence on x 
in the Lagrangian, although these symmetries are spontaneously broken by the kink. Starting from one 
static kink located at the origin, we obtain the complete family of kink/anti-kink solutions: 

(t)K{x) , x = {-lf—. '^^~^ ^ , where a = 0, 1 . (10) 

V 1 — v^ 

Although the coordinate (10) x should be entered in all the formulas in the sequel regarding operators, 
superpotentials, kinks, etcetera, we shall always merely write x (meaning that we set the kinks at their 
center of mass and identify kinks with anti-kinks) in order to alleviate the notation. 
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2.2 Kink stability and super symmetric quantum mechanics 



The field small fluctuations over the kink background 6f{t,x) = (p{t,x) — (pxix) are still solutions of the 
second-order field equations (2) if the linearized PDE 



6f{t,x) = 0{6') , a = ^-^ . (11) 



is satisfied. In this (static kink) case, the general solution can also be obtained by means of separation 
of variables, 6f{t,x) = eyip[iujt]f^{x), but one needs to solve the spectral problem KN-fu}[x) = uj'^fui{x) 
of the Schrodinger operator 



Kn- = -:n?+vl+ Vn-{x) where Vn-{x) = -jT?[(t>K{x)] - vl (12) 



instead of relying on the eigenvalues and eigenfunctions of a Helmholtz operator as in the vacuum case. 
Here the natural number € N refers to the number of bound states of the V/v- (x) potential well. We 
shall restrict ourselves to models such that the function Vn^{x) tends to zero at both x = ±00 and, 
moreover, we also require: 

Thus, the particle masses emerging around the two vacua connected by the kink have the same value: 
= u^. Bearing this in mind, the asymptotic behavior of (12) is equivalent to (IS)"'. 
The stability of the solution cpx is guaranteed when the eigenvalues of the operator (12) are non- 
negative, because in this case the fluctuations Sf{t, x) remain bounded in the temporal evolution. In 
general, the spectrum of the operator (12) is of the form 

Spec{KN^) = Wl}n=0,l,...,l-1 U {Ujf}s^ U + vl}km , UJ^ < Wn+l , (14) 

where a;^ are the bound state eigenvalues, / + 1 is the number of eigenfunctions in the discrete spectrum, 
and Lo'^{k) = k"^ + v'^ correspond to the continuous spectrum eigenvalues of the Ktv- operator. According 
to the one-dimensional Levinson theorem, the subscript si in (14) to be entered in the DHN formula, 
distinguishes whether the highest eigenvalue in the discrete spectrum, cof, comes from a bound or a half- 
bound state. If ujf < this eigenvalue corresponds to a (normalizable) bound state, = / + 1, and we 
shall choose s/ = 1 in the DHN formula, whereas if uif = v'^ there is a (non-normalizable but bounded 
at infinity) half-bound state just at the threshold of the continuous spectrum; thus N = I, and the DHN 
choice will be s/ = |. 

To show the kink stability we differentiate with respect to x the ODE equation (7) evaluated at the 
kink solution: 



d d^(f)K d 5U,^ , _ ( d^ 



^)] ^ - ( + ^-^IMx)] ) ^ = . (15) 



dx dx'^ dx 6(1)^^^^^ " V dx"^ ' 5(j)^^^^^^ ") dx 

The wave function 



Mx) = ^{x) (16) 
dx 



^There are interesting models where the vacua connected by the kink are not equivalent, giving rise to particles with 
different masses, e.g., the Khare-Lohe cjP model. One understands that the moduli space of vacua is formed by several 
orbits of the symmetry group. The quantum mechanical spectral problem posed by the kink Hessian in this situation is non 
unitary. 
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is accordingly an cigcnfunction of the K^^ operator with eigenvalue zero. For topological kinks, mono- 
tonic functions of x, their derivatives have no nodes and the translational mode is the ground state 
of the second-order kink fluctuation operator, i.e., ujq = 0, and topological kinks are always stable. Non- 
topological kinks, however, give rise to translational modes with one node and cannot be the ground 
state. There is an cigcnstatc of negative energy and non-topological kinks are unstable. Thus, the one- 
loop mass shift of non-topological kinks picks an imaginary part since it corresponds to resonances in 
the quantum domain, whereas stable topological kinks are bona fide quantum eigenstates of the field 
theoretical Hamiltonian. 

The existence of the zero mode ground state is the crucial fact for building the structure of super- 
symmetric quantum mechanics in the K^- spectral problem. The Kj^-. operator can be written as the 
product of two first-order operators in the form 

K]s[-=A^j^An with An = -^+Wn{x) and A^j^ = --^ + Wn{x) , 

where the function Wn{x) is minus the derivative of the superpotential 'Ej>^{x) = log V^oCa^); the logarithm 
of the ground state wave function: 

^.W = -H'.W = -||| . (17) 

Explicitly, the second-order kink fluctuation if at- operator is defined in terms of the Wn{x) function: 

KM- = -^ + WN{xf -W'n{x) , WN{xf-W'i,{x) = vl+VN-{x) . (18) 

Quite generically the asymptotic behavior of the derivative of the superpotential at very large distances 
from the origin is : Wn{x) ~a;_>.±oo (~l)"^jv tanh[i;jyx] where a = 0, 1. Thus, 

lirn Wn{x) = ±{-ITv^ , lim W'j^{x) = . (19) 

x—^±oo X— >±oo 

There is a partner Schrodinger operator that is also factorized in terms of the first-order An and ^jy 
operators but in the reverse order: 

Kn+ = AnA^n = + WN{xf + W'n{x) = -^+vl + Vn+{x) . (20) 

All this calls for the construction of the following M = 2 supcrsymmetric quantum mechanical system: 
there are two supercharges 



\ _ / A 



that close the supersymmetry algebra 



Q% = iQ^Nf = , Kn = QnQ^n + QnQn ■ (22) 



The Hamiltonian 
Kn = 



Kn- \^ ~£2+WN{xf-W^{x) \ 

Kn+ J [ -£^ + WNixf + WUx) J 

is supersymmetric because it commutes with the supercharges, KnQn ~ QnKn = K^Q^n ~ Q^n^n = 0, 
which are consequently constants of motion of the system. In sub-Section 2.4 we shall describe the very 
well known properties of the spectrum of supersymmetric Hamiltonians that will help us in our main 
task of computing one-loop kink mass shifts. 
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2.3 SUSY quantum mechanics, the BPS equation, and reconstruction of scalar field 
models supporting kinks 




We now shall explain how the supersymmetric quantum mechanical structure of the topological kink 
stability problem is indeed deeply encrypted in the BPS or first-order ODE (9)^. Plugging the field 
configurations written as the kink plus a small fluctuation (p{x) = (pxix) + 5(l){x) into the BPS equations 
we find new solutions if and only if the following linear ODE holds: 



x) = . (23) 



Because the topological kink solves equation (7) we know that U'{(j)K) = V'o the equation (9) means 
that ^J2U{^^k) = i'o- The differential operators in the equation (23) are no more than the factorization 
operators An and ^jy. Therefore, the factorization property of the second-order kink fluctuation operator 
is closely related to the existence of first-order equations in the scalar field theory model. 

This, at first sight, hidden property of BPS kinks is in the backyard of the mechanism essentially 
discussed in reference [13] and used explicitly in papers [36, 37]. Starting from a given superpotential 
function Eiy{x) the field theory model can be reconstructed by following the subsequent steps^: (a) From 
the equation (17) we write the zero mode in terms of the superpotential: 

V;o(x) = e=^(^') = e- -^^ '^^^^^^^ . (24) 

(b) From the zero mode (16) the kink solution is obtained through indefinite spatial integration: 

ct>K{x)= r Mi)di . (25) 



(c) The final step is the most delicate one. The first-order equation (9) allows us to write the function 
U{cj)) characterizing the field theory model 

U{<t>) = \i^l{rK\(t>)) , (26) 

from the inverse of the bijective function (j) = tpxix), e.g., x = log tan ^ for the sG-kink, or x = arctanhc/) 
for the <f>^ kink. Formula (26) defines U{(j)) in the range of the kink cI)k{x), i.e., in the interval (l){t,x) G 
[(j)y\ (py'^^'^]. Some type of C^(M)-continuity is needed to extend ?/(</>) to the rest of the field space. We 
shall make use of this procedure in the section 3 in order to construct some field theoretical models 
starting from a given kink zero mode. 



2.4 Kink mass quantum correction from the SUSY partner operator 

After identifying a classical kink solution we are interested in analyzing its quantum properties. Pursuing 
this goal Dashen, Hasslacher and Neveu succeeded in computing the one-loop quantum correction to the 
classical mass of the kink and soliton found in the (p"^ and sine-Gordon models [1, 11, 12]. After a 
zero-point renormalization and a subsequent mass renormalization, these authors obtained what is now 
known as the DHN formula, which has been generalized in ['2 !] for operators having a non-transparent 
scattering spectrum. The generalized DHN formula applicable both to transparent and non-transparent 

*This structure has been derived in the previous sub-Section from the existence of the translational zero mode. 
^Another interesting procedure of reconstruction based in the inverse scattering method has been proposed on [38]. 
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Hamiltonians reads: 



d6(^~Hk) 1{Vn-{-)) 



dk 



N 



N 

(27) 

where are the bound state eigenvalues of the second-order kink fluctuation Kn- operator, si is equal 
to ^ if the spectrum of the Kn- operator involves one half-bound eigenfunction and it is one otherwise^. 
Moreover, S^^~\k) denotes the phase shifts in the scattering waves of the K]\f- operator and {V]\i-{x)) 
stands for the area enclosed by the function Vn-{x) and the real line: (VN-ix)) = V]y-{x)dx. We 
recall that Vj\[-{x) is the potential well of Kj\[- minus the non-zero threshold v'j^. It is worthwhile at this 
point to explain the physical origin of the different terms in formula (27): 

1. The first two terms encode the contribution of the singlet (bound and half-bound) states to the 
kink fluctuation energy. The third term corresponds to the subtraction of the singlet state of the 
vacuum fluctuation spectrum. 

2. The positive term in the integral over the scattering modes collects all the kink minus vacuum 
continuous spectrum fluctuations to the energy. 

3. The finite fourth term comes from using an intermediate regularization setting the same number 
of kinks as vacuum fluctuation modes, rather than a cutoff in the kink and vacuum energies. All 
these terms assemble together make what one might call kink Casimir energy. 

4. The remaining, negative, term in the integral is due to the contribution to the fluctuation energy 
of the one-loop renormalization mass counterterm. 

Nevertheless, the use of expression (27) brings with it a handicap: complete analytical knowledge 
of the spectral data concerning the second-order kink fluctuation Kj\[_ operator is a rare occurrence in 
generic (l-l-l)-dimensional scalar field models. This makes the computation of the kink mass shift seldom 
feasible. The goal of this paper is to take advantage of the techniques available in ID supersymmetric 
quantum mechanics in order to alleviate the dependence on the spectral information of the DHN formula 
(27). 

The aim in this section is to rewrite the generalized DHN formula (27) in terms of the spectral data of 
the supersymmetric partner operator to K^^. The supersymmetric algebra shows that the spectral 
data concerning the Kjsr- operator, defined in (12) and (18), and the data of its supersymmetric partner 
Kn+, defined in (20), are closely related [8]. First, bearing in mind that V]\f+{x) = Wiyix)"^ + W'j^{x) — vj^ 
we find that 

lim Vn+{x) = , 

>±oo 

showing that the K^- and Kn+ operators share the threshold vf^ of the continuous spectrum, see Figure 
1. Second, if the spectrum of the -f^Ar- operator is summarized in the form (14) then the spectrum of 
ifjv+ reads: 

Speci^Tn- = {ujl}n=i,...j-i U {ujf}si U {k"^ + v%}j, 



Both spectra look identical but there are two important differences: (1) In Speci^Tn- there is no zero 
mode. (2) Although the continuous spectrum eigenvalues are the same the spectral densities are not. 



®We recall that a half-bound state is a singlet eigenstate of energy eigenvalue equal to the energy of the threshold 
continuous spectrum. The reason for the | factor (and the name) is the one-dimensional Levinson theorem. 
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t Vif+Vjf+ix) 







Figure 1: Graphical representation of the spectra of the kink fluctuation operator Kjy- and its SUSY partner 
operator K^j^ (red horizontal lines) and the potential well arising in these operators (solid blue lines). 



Third, the phase shifts emerging in the scattering eigenfunctions of these two K]\f^ and Xat- operators 
are also related [ ] . Given the asymptotic behavior of the right-going scattering waves of Xat- and i^Ar+ 

V'(^+)(A;,x^ -oo) ^e^^^ + pS^^e-*'^^ , ^^^+\k,x ^ oo) ^ a'}:^_^e''''' , 
where A:^ = cj^ — vf^, the relations between the Kj\f± reflection and transmission scattering amplitudes 



are: 



WN-+ik (^) (^) WN+-ik (^) 

PN- = ^^^^—f^PN+ ' "^-=W^vr^^^+ ' -hereTF;v± = ^lim^W^^(x) . 

Analogously, for an incident wave coming from the right side we have the relations: 

(^) _ Wn+ - ik (4_) _ Wn~ + ik (^) 

" H^jv+ + ik ' " Wn+ + ik ''^+ ■ 

Without loss of generality we choose W]\f+ = = —Wn-- In this situation cr^^_^ = f^"* = crN+ and 
= <7m1! = <7Af-- The identities immediately above reduce to 



(^) _ V^-ik (^) (^) _ V^-ik _ v^-ik 

Pn~-^;^Pn-, , Pn--^;^Pn+ > - • (^8) 

The unitary scattering matrix Si\f- encoding the scattering data of the operator /Cat- can also be defined 
in terms of the reflection and transmission scattering amplitudes of the supersymmetric partner operator 



S 



Pn^ J -n^ 



The phase shifts 6^ ^ are the arguments of the eigenvalues A^i^^ = e^**^* of the matrix Sn- . Thus, 
the relation between the phase shifts of the supersymmetric K^^ / Ki\f^ operator pair follows: 



lis^'^'^ , / (-5.) u„ — ik 



, (^) (^) 
aN+T y Pn+Pn+ 



.!^^24--) . (29) 
v^+ik ^ ^ 



Defining a relative phase shift by the formula e ° = i.e.. 



2-y /j; 

^o(fe) = arctan „ ^ „ , (30) 
— vt 
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expression (29) becomes 

Therefore, the connection between the Kn± total phase shifts 6^^^\k) = 5^_^^\k) + 6[!^^\k) is 

6^^-\k) = 6^^+\k) + eo{k) + TT . 

The generalized DHN formula (27) requires the derivative of this phase shifts with respect to the mo- 
mentum k, which is the spectral density of the kink fluctuation minus the spectral density of the vacuum 
fluctuation, 

dd^^-^k) _ d5^^+\k) _ 2vN 

dk dk k"^ + v'j^ ' 

where we have used (30). Plugging this result into (27) we obtain a new version of the generalized DNH 
formula 



k'^ + V 



(31) 

Fortunately the 5^^^^^ (k) phase shifts are related by supersymmetry in such a way that the modification of 
the DHN formula merely requires a shift in the mass renormalization contribution to the kink fluctuation 
energy. The new formula depends on the spectral data of the supersymmetric partner operator i^Ar+ 
instead of K^- itself. Sometimes the analysis of the spectrum of is simpler than that of Xat- and 
the 5^^~^\k) phase shifts are analytically accessible, whereas the 5^^^\k) either are not accessible or 
are much more involved. In this situation formula (31) is more manageable than the original generalized 
DHN formula (27). 



2.4.1 One- loop mass correction from the sine-Gordon soliton fluctuations 

We shall illustrate the use of (31) by means of its application to the computation of the one- loop quantum 
correction of the sine-Gordon soliton mass. In this model we need only the A'^ = 1 version of the general 
formalism. Use of formula (31) renders this calculation almost trivial. The sine-Gordon potential is 

U{(j)) = 1- cos (j) . (32) 

The degenerate absolute minima of f/(<^) form the discrete set A4 = {(j)^^^ = 2iTj}, j € Z, and hence the 
masses of the fundamental quanta are: vf = |^ = l,Vj. Besides these static homogeneous solutions 
the model admits the kink extended solutions 

(psix) = ±4arctane^' + 27rj , 

which interpolate respectively between the (p^^^ = 2-Kj and (f)^^^^'> = 27r(j it 1) minima. The second-order 
sG-soliton fluctuation operator (12) is 

Ki_ = -— + l-2sech^x , (33) 

which fixes Vi-{x) = — 2sech^3; and therefore {Vi-{x)) = —4. Following the scheme described in the 
section 2.1 we find from (16) that the zero mode of Ki- is: ipo{x) = 2 sech x. Also this Hamiltonian admits 
a half-bound state: V'i/2(^) = tanhx of energy one. Therefore, the superpotential and its derivative are: 
Hi(x) = — logcoshx, Wi{x) = tanhx, as one may see from (17). The = ^ +tanhx first-order 



11 



operator permits the factorization of (33) in the form Ki- = A\Ai and unveils its supersymmetric 
partner (20) : 

^'- = -^"^! = -|^ + ' ■ 

Obviously Vi^{x) = and the phase shift is null: 5^^^\k) = 0. The spectrum of the Ki^ operator 
includes the continuous spectrum of plane waves and energy uP{k) = k"^ + 1 plus one singlet (half-bound) 
state, a constant function, of energy one. Therefore, there is only one bound state in the spectrum of 
the SUSY pair: the Ki^ zero mode. The one-loop sG-soliton mass quantum correction is computed by 
applying the formula (31) to find 



K) 



dk 



-2 + 2 



1 

TT 



which recovers the classical result of Dashen, Hasslacher and Neveu [1, 11] in an elementary way. 



2.5 Kink mass quantum corrections from SUSY quantum mechanical hierarchies 

In the previous section we profited from the existence of a supersymmetric partner K]sr^ to the second- 
order kink fluctuation operator K^- in order to obtain a new DHN formula (31), which depends only 
on the spectral data of K^^. We shall continue exploiting this strategy by building a SUSY quantum 
mechanical hierarchy that will allow us to reduce the kink mass computation in terms of the spectral 
data of simpler Hamiltonians in many more models. Starting from the (as yet unspecified but factorized) 
operator Kn^ = ^jv^jv we define a new operator as 

d^ 

K(N-i)- = Kn+ -un = + v%-u% + Vn+{x) = + W(^7v-i) + ^iN-iyix) , 

where ^^^^_;^-) = ^^at ~ ""at ^-^cl V(Ar_i)_(a;) = Vn+{x), as the starting operator of a new stage. If the 
-Rr(Ar-i)- operator admits the factorization 

K(^N-i)- = ^(Ar-i)^(W-i) with = ^ + W(^N-i)ix) and A\^_^^ = + ^{N-i)^^) . 

the strategy explained in the previous section can be applied iteratively. This happens with the choice 
uff = Lof, i.e., we subtract the energy of the first excited level of K^^ identical to the ground state energy 
of K]\f^. In terms of the new function W(^]\f_i-){x) the K(Af-i)- operator reads: 

KiN-iy = + W^N-i){x)^ - ^(V-i)(^) • 

Taking into account that the definition of ir(Af-i)- only involves a translation of the Ki\i^ origin of 
energies, the spectra of i^(Ar_i)_ and K]\f^ are related trivially. Denoting 

SpecK(Ar_i)_ = {f^}„=l,...,i-l U {ui}s, U {A:^ + W^Af-l)}fc6]R 

the relation i'^ = uj'^ — = u'^ — (v^ — w^^ -^p shows that z/f = and the kernel of ii'(Af-i)- exhibits 
a zero mode. The phase shifts associated with the -ftr(Af-i)- and K^-^. operators obviously coincide: 
<5(^-i)-(fc) = <5^+(fc). 

The following step is to construct the supersymmetric partner operator to -ftr(Af-i)- 

^2 ^2 
K^N-l)+ = ^(7V-l)^(^_i) = + W(^M-l){xf + VF(V_i)(x) = + vf^_^^ + V^N-l) + {x) , 



12 




Figure 2: Graphical representation of the energy spectra of the SUSY-partner operator hierarchy K^zp, K(^n-i)^, 
(red horizontal lines) and the potential wells characterizing these operators (solid blue lines). 



where the notation V(jv-i)+(2^) = W{N-i){x)'^ + ^(V-i)(^) ""^^^at-i) evident see Figure 2. Again, there 
exists a connection between the spectral data of the and Ki\i^ operators: 

Speci^(Af_l)+ = {iyl}n=2,...,l-l U {l^f}s, U {k^ + t;^Ar_i)} , 

note the absence of the zero mode, and the 6^^^~^^'^\k) phase shift derivatives satisfy: 



d5{(^-i)+)(fc) d5((^-i)")(A:) 



2v^ 



dk 



dk 

d5^^~\k) 
dk 



+ 



(7V-1) _ dS^^'+Xk) 



2v 



A;2 + ' 



+ 



2vn 
k^ + v% 



'(N-l) 

2v 



dk 



+ 



(7V-1) 



/t2 + 



'{N-l) 



(N~l) 



Plugging this last result into formula (31) the generalized DHN formula is promoted to this level of the 
hierarchy: 

i-i 



AE{(I)_ 



•K 



1 



n=l 



+ 



27r 



dk 



d5({^-l)+)(A;) 



dk 



k^ 



UVn-{x)) + 2vn '^^(N-l) 



k^ + 



'N 



+ '^In-1) 



(34) 



Note that this depends on the spectral information of the iir(7v-i)+ operator. 

Formula (34) can be further generalized by repeating these steps in a recursive way. By induction 
it is easy to check that the phase shift derivative arising in the -ftr(7v-j)+ operator, the j-th level of the 
supersymmetric hierarchy, satisfies the condition 



d5((A^-i)+)(fc) _ d(^((^-i)-)(A;) 2t.(^-i) _ d5^^-\k) 



dk 



dk 



+ 



fc2 + 



{N-j) 



dk 



j=0 



2v, 



(N-i) 



'(N-i) 



where f .) 



'N 



l-l 



Lof. In this level of the hierarchy the generalized DHN formula 



1 v-^ 1 
9 2^ ^"+9 



+ 



n=l 
1 

2^ 



sm - X + ^ (^iV-(^)) + 



dk 



d5((^-j) + )(A;) 



dk 



A;2 



\{Vn-{x)) + 2vn 



J 2v 



'N 



'N 



E 

i=l 



(N-i) 



k^ + 



'N 



k'^ + "^fN-i) 



(35) 



expresses the kink mass quantum shift in terms of the spectral data of the j'-th supersymmetric operator 
with trivial kernel. 
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2.5.1 One-loop quantum correction to the (/!)^-kink mass 



The 4>^ model is characterized by the potential energy density (/{(p) = \{(t)'^ — 1)^ in the action functional 
(1). The set of vacua is M = {—1,1}. The fundamental quanta emerging from any of these two vacua 
have masses equal to u| = 4. The kink solutions 

4'k{x) = tanhx 

connect the two vacua of the model. In this case N = 2 because the 0^ kink second-order fluctuation 
operator presents two bound states: 

K2- = —-^ +4 — 6 sech^ X 

The potential well V2- (x) = —2(2 + 1) sech^ x is the second in the hierarchy of Poschl- Teller reflectionless 
potentials. Moreover, {V2-{x)) = —12. From (16) and (17) the zero mode is identified as tpoi^) = 
sech^ X, whereas the superpotential, its derivative and the factorization operator are respectively: ^2(3;) = 
— 21ogcoshx, W2{x) = 2tanhx, and ^2 = ^ + 2tanh a;. The supersymmetric partner operator to 
K2- = A2A2 is also reflectionless and Posch- Teller, such that 

K2+ = A2AI = -— + 4 - 2 sech^ X = A\Ai + 3 = Ki_ + 3 

defines the next level in the supersymmetric hierarchy. The Ki- operator in the formula above in turn 
factorizes in the form = where Ai is the first-order operator arising for sine-Gordon kinks. 

We have implicitly used the fact that o;^ = 3 is the first excited eigenvalue of K2- and the ground state 
energy of i^2+- The supersymmetric partner operator to Ki- is, as we know from the sine-Gordon kink 
analysis, 

dx"^ 

which fixes vf = 1. The phase shift associated with this free particle operator is null: 5^^^\k) = 0. 
Again, the SUSY quantum mechanics hierarchy allows us to apply formula (34): 



Ki+ = AiA\ = -— + l 



AEUk) I ^ 12 1 



VS - — + — dk 



hY^ 2 All 2tt 



00 







-6 + 4 2VA;2 + 4 



2 vr ^ 2^/3 vr 



The one-loop kink mass quantum correction obtained this way is exactly the classical result of Dashen, 
Hasslacher and Neveu [I, 11, 12]. 

2.5.2 One-loop kink mass correction in the parent potential models 

In references [6, 13] the authors introduced and discussed a hierarchy of interesting (l+l)-dimensional 
field theory models (which were referred to as parent models) whose action functional (1) is specified by 
the potential terms 



2r[i + f]2 



^^^^W = ^^^LVJ a-^'MH;if])'' • (36) 



vrr[f]2 



Here, A'^ G N is a natural number and / ^ [x; a, h] is the inverse of the regularized incomplete beta function 
(or regularized beta function) I{x] a, h). The authors also identified the kink solutions in these models 

(jl-p (x) = Sign(x) /[tanh^ x; i , f ] , (37) 
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all of them connecting the vacuum points (j)v = ±1- The second-order differential operators governing 
these kink fluctuations are precisely the reflectionless Poschl- Teller Schrodinger operators shifted in such 
a way that the ground state energy is zero: 

(P 

Kp^^ = -— + N'^ -N{N + l)sech^x . (38) 

In fact the motivation in these works is the identification of the field theory models whose second- 
order kink fluctuation operator is the Schrodinger operator (38). These operators have been studied 
in depth in the literature. The discrete spectrum of K^^^ is formed by the discrete set of eigenvalues 

= n{2N — n), n = 0, 1, . . . , — 1, together with a half-bound state, a;|^ = N"^, at the threshold 
of the continuous spectrum, see e.g. [39, 40]. The phase shifts, however, are defined in terms of very 
complicated hypergeometric functions. This is the point where we shall take advantage of (35). This 
formula computes the kink mass quantum correction in terms of the phase shift of any operator belonging 
to the supersymmetric hierarchy that, in these models, ends in the free Helmholtz operator where the 
phase shifts vanish. The application of formula (35) to obtain the mass shifts due to the kink (37) 
fluctuations will be thus straightforward. 

Prom (38) we read = A^^ and Vn-{x) = —N{N + l)sech^x, such that the area enclosed by this 
function is {Vn-{x)) = —2N{N + 1). The superpotential and its derivative, in turn, are respectively 

Siv(a^) = —iV log cosh X and Wn{x) = A^tanhx. K^- = A]^An factorizes in term of the first-order 

A 

dx 



operator: A]\[ = 4^ + A^tanhx. The supersymmetric partner operator is thus: 



Kn+ = -j^ + {N-l)'^ -{N- l)A^sech2 x + 2N - 1 = + 2N - 1 (39) 



We also recognize in formula (39) the operator arising in the next level of the supersymmetric 

,2 



hierarchy and identify v?J^_^^ = {N — 1)^ as well as Vn+{x) = V(Ar_i)_(x) = — (A^ — l)A^sech^x. The 



first excited bound state eigenvalue of K^- is iof = vf^— = N'^ — {N — 1)^ = 2N — 1, as requested. 

Applying this process iteratively, the hierarchy of the supersymmetric pairs of operators is constructed. 
The superpotentials and their derivatives in the j-th iteration are: H(;Y_j)(x) = — (A^ — j) log cosh x, 
VF(Ar_j)(x) = (A^ — j) tanhx. In the next array we collect the pertinent information 

KiN^,)- = + - •?')(^ - + l)sech2x , vf^_^^ = {N- jf 

K^N^j)+ = -^ + {N-j-lf-{N-j-l){N-j)sech^x + 2{N-j)-l = 
= i^(jv-j-i)_ + 2(Af - j) - 1 

In the last iteration, j = A^ — 1, we find that = —-^ + 1 and therefore 6^^~^\k) = 0. Plugging the 
j = N — 1 information into equation (35) we obtain 



dk 



N{N-1) 2t;^^_^.yfc2 + Ar2 



Because ^jlf/ (A^ - j) = 2N{N - 1) we write the previous formula in the form 



1 r^^r ^ N(N + 1) 1 , v-^ 

n=l j=l 



Vk^ + m 



\ - - I X; ^/JW^) arccos ^ 

n=l j=l 
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1 1 ^""^ N - ■ 

= -—N{N + 1) + - V Vj(2iV - j) arcsin — — ^ . 
Ztt tt — ' iV 

i=i 

The one-loop kink mass correction in the parent models (36) is finally 

^^f^^^fJ. = -±N{N + 1) + -Y ^/W^ arcsin , (40) 
27r vr ^ AT ' ^ ^ 

" r=l 

where we have redefined r = — j in the summation index. The same result was found by Boya and 
Casahorran in [14] by applying the procedure of Cahill, Comtet and Glauber [ i ], which gives the one- 
loop kink mass correction in terms of the bound state eigenvalues of the K^- operator, see [14]. We 
stress that this last procedure is valid only if the potential well is reflectionless. 



3 Exotic scalar field theory models with exotic kinks 

In this section we shall discuss other scalar field theory models, where the issue of the computation of the 
one-loop correction to the kink masses is an analytical nightmare if one uses the standard approaches to 
the problem. We shall show, however, that in these exotic cases the supersymmetric quantum mechanical 
techniques works well in computations of the kink mass shifts. 



3.1 Kink mass quantum shift in the U{(p) = |</)^ cos^ log |0| model 

The action (1) in this model is specified by the potential energy density 

[/(</)) = ^Os2 log |</.| . 

This function is depicted in Figure 3. Notice that U{(j)) as a function of (/) oscillates indefinitely around 
the origin. We check, however, that U{(j)) G C°°(M - {0}), although only U{(p) G C^(M) in the whole real 
line. This model was addressed in reference [• . :], where the authors succeeded in identifying the kink 
type solutions. 





Figure 3: Re-scaled graphics of U{(j)) (left) and plot of the kink 



^Ki connecting the vacua 



and 



(right). 



The set of vacua M of the model is the set of zeroes of U{(l)): (jy'y 







log|(/)^^'^ 



V 



7r(j + i),iGZ. 



Therefore, the set M is the union of two subsets separated by the origin which is a singular point: 
Thus, there is an infinite number of vacua in this model: 



1 I^^UI ff.TJ 

2(l + cos(log<A^)-sin(log<^2)(3 + |</>|^)) , |^|^(,) 



41 



1 , Vj 
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all of them being equivalent in the sense that the masses of the fundamental quanta emerging from any 
1. The distance in field space between consecutive vacua, however, differs with j: 



vacuum are: v 



V± 



62 



+7rj ( -K 



(e--l). 



The kink solutions of the first-order equation (9) are 



4\{x) 



•gd{x)+TTj 



arctan tanh ^ +7rj 



„2 arctan tanh ; 



-irj 



(41) 



which respectively interpolate between the (/>y^ and (kv^^ vacua In Figure 3, a particular kink 
connecting the (l^y^ and the (/)y^ vacua is depicted. 

In this section we shall compute the one- loop correction to the classical mass E[(l)"^_^ = fe^'^-' coshvr 
of these kink solutions (41). The second-order operator governing these kink fluctuations is 



+ 1 — sech X — 3 sech x tanh x 



-2. 



This sets Vi-{x) = — sech^x — 3sechxtanhx, which is the Scarf II potential, and therefore (yi-{x)) 
The zero modes are identified from condition (16) 

[x) = sech X a^-^taiita^h f +7ri 

and by using (17) the derivative of the superpotential is Tyi(x) = tanhx — sechx. This is a remarkable 
fact: despite being of different "length "the kinks labeled by different j give rise to identical wells*^. The 
second-order kink fluctuation operator Ki^ factorizes in the form 



Ki_ = A\Ai with Ai 



dx 



+ tanh X — sech x and A 



A. 
dx 



+ tanh X — sech x 



and consequently the supersymmetric partner operator to Ki^ reads: 

d^ 



K 



1+ 



AiA\ 



dx"^ 



+ 1 + sech X — sech x tanh x 



Although the Ki- and operators look similar, the analytical expressions of the eigenfunctions in 
the continuous spectrum of are much simpler than those associated with the operator. This 
allows us to extract the spectral data, the bound state energies and the phase shifts of this operator, 
in an easier form and reduce the calculations required to obtain the kink mass quantum shift. Figure 4 
shows the potential wells arising respectively in the i^i- and Ki+ operators. 




n-+ii_(x) 




Pa 


i-f + I'l+lx) 




.V 


-10 -s 

-1 


i 10 



Figure 4: Potential wells of the (left) and K\j^ operators (right). 

gd(a;) is the standard notation for the Gudermannian. 
®The connection of the (inverse) Gudermannian to the Mercator projection is well known. Lengths are distorted in such 
a way that different parallels look similar: the parallels closer to the Poles grow longer. The kinks in this model exhibit a 
similar phenomenon but the rescaling is exponential. 
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The eigenfunctions in the continuous spectrum of Ki+ are 



, , 1 1 + zsinhx" 

— ik, ik 



2 ' 2 



1 1 3 1 + i sinh x 

+ i — ik, - + i + ik, - +i. 



-2 2 ' 2 ' 2 

where = + The asymptotic behavior of these two hnearly independent eigenfunctions for x — oo 
is 

2Vvrr[i - i + ifc] 2Vvrr[i - i - iA:] 

w„fa;) ""^^ ^-f(fc-i) + ^fc]^[| + z] ^g(fc+i) - ^fc]^[| + z] ^.^^ 

^ 2V^A;r[^ + i + zA;] 2^/^kT[l + i - ik] 

whereas for x — >■ +oo it becomes 

^zis? gg(fc-i) + zA;]r[| + z] ^^^^ ^-f(fc+i) - + z] 

2V^A;r[i + i + iA;] 2V^A:r[i + i - ifc] 

A Hnear combination of the continuous functions V'i(^) 'tp2{x) describes an incident wave coming 
from the left side, behaving asymptotically as: 

^ik. , 7re-fe-r[l + ik] (sech (fc + l)7r - e^^ sech (1 - k)n) ^ ^ -oo 

{1 + e'^'')r[l - ik]r[l - i + ik]r[l + i + ik] ^ x oo 

7re^r[^ + ik] (sech {k + 1)tt + sech (1 - k)7r) .^^ 
(1 + e2^)r[i - iA;]r[i - i + iA;]r[i + i + ifc] ^ x ^ +oo 

This result determines the reflection and transmission scattering amplitudes of this incident wave: 

(^) (t.] _ ^e-'=''r[^ + ik] (sech {k + l)7r - e^'^ sech (1 - k)Tr) 

^ ^ ~ (1 + e2'^)r[i - ik]T[l - i + zA;]r[i + i + ik] ' 

(_^) 7reT[i + ik] (sech (A + l)7r + sech (1 - k)^) 

^ ' ~ (1 + e^^)T[l - ik]r[l - i + ik]r[l + i + ik] 

In a similar way, the asymptotic behavior of an incident wave coming from the right side is 

Tre'Til + ik] (sech {k + l)7r + sech (1 - k)7r) .-^^ 



(1 + e^^)r[l - ik]T[l - i + ik]T[l + i + ik] 



X — >■ — oo 



ire'^-nl + ik] (e^^sech {k + l)7r - sech (1 - k)7r) x ^ +oo 

{I + e'^^)r[l - ik]r[l - i + ik]T[l + i + ik] ^ X oo , 

such that the reflection and transmission scattering amplitudes are: 

(^) _ Tre''''T[l + ik] (e^'^sech (A; + l)7r - sech (1 - k)Tr) 

" (1 + e2-)r[i - ik]r[l - i + ik]T[l + i + ik] ' 

(^) _ 7reT[^ +i/c] (sech(A + l)7r + sech(l - A;)7r) 

"^^^ ^ ^ ~ (1 + e2'^)r[i - ifc]r[i - i + iA]r[i + i + iA] 
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Note that a^j^\k) = a[^\k) = ai^{k) and p^^\k) = —p^^\k). The unitary scattering matrix follows: 

5*1+ = 




The phase shifts are the arguments of the modulus-one eigenvalues of this matrix 

2i5i^+) , • (^) + ik] 

— i^]r[2 — i + ik]T[2 + i + i/c](sinh7r it icosh/cvr) 

The derivative of the total phase shift 6^^^\k) = S^^\k) + 6^^^\k) with respect to the momentum k 
becomes 

^^^^^ = Re[^{l-ik)+^{l+ik)-^{l-i + ik)-iP{^ + i + ik)] , (42) 

where ^(z) is the digamma function. We stress here that the calculation of the spectral density coming 
from the scattering waves is much easier to achieve than if we had handled the corresponding 
computation for the Ki- operator. Note that Ki- has only one bound state, the zero mode, because 
there are no bound states in the spectrum of Plugging (42) into formula (34), the one- loop correction 
to the (41) kink masses can be figured out: 



AEi(l)K) 



2 



1 1.x 1 /■ 

4 + 4^-<-^)+2^X 



OO 

dk 



dk Vk^Tl 



This expression requires the integration of a complicated function, which we estimate numerically to find: 
dk['^^^-^Vl^Tl - ~ -1.1823055. Therefore, the final response to the mass quantum shift 

of all the kinks in this model is 

^t^^ ^ -0.597325 . 
3.2 A model reconstructed from a zero mode 

The procedure proposed in reference [ ] shows how to construct a field theory model from a specific 
function selected as the zero-mode of a conjectural second-order kink fluctuation operator. This procedure 
has been used by Kumar in order to generate models with very general soliton profiles, see [36, 37]. The 
steps to be followed in order to implement this idea are described at the end of section 2.1. Here, we 
start from a particular selection of the zero- mode (see Figure 5, left): 

^o(x) = e'g'^^^) sechix = 6-2'^^'=*^°*'^"^ | ge^hlx , 
This function is the ground state of the following second-order fluctuation operator 

d^ I 1 

Ki- = — — T -|- - -|- T sech^j; -|- 2 sech x tanh x . (43) 
<ix2 4 4 

The kink proflle is identified through the integration of the zero mode, (25): 

Here, the change of variable, z = C(x), where ^ : M — )■ (— f , f ), x i— ?■ ^{x) = 2 arctan tanh |, has been 
performed for the sake of simplicity. The inverse of this diffeomorphism is x = ^"^(z) = 2 arctanh tan | 
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Figure 5: Zero mode (left), potential well (middle) and kink (right) encompassed in the second-order kink fluctu- 
ation operator Ki^. 



and, in particular, dx = sec zdz. We also define the transcendent function T : [— |, |] — )■ M as T(z) = 
e~^\/ sec tdt such that the kink solution 

(j)Kix) = T(2arctantanhf) (44) 
is expressed in terms of the hypergeometric function in the form 

c^k{x) = Re {(I - |)e-2(i+^)-t-t-h| ^/^^ (^-1 + [-{ + ^,1, ^ + • 
The classical kink energy, however, is easy to compute: 

/oo roo 
^o(x) dx = exp[— 2gd(x)]sechj; = sinhvr . 
-oo J — oo 

Finally, the potential energy density U{(p) that characterizes the scalar field theory model is determined 
from (26) in the kink profile range: 

Uic^) = ^e^^''^^^cos[r-\cl>)] . 

The second-order kink fluctuation operator (43) sets the values vf = ^ and Vi- (x) = ^ sech^x -|- 
2sechxtanhx. Hence {Vi-{x)) = ^ holds and the function VFi(x) allowing the factorization of the 
Ki- operator is VFi(x) = —-£:\ogipo{x) = |tanhx -|-sechx. In order to estimate the mass quantum 
shift we need to extract the spectral information of the operator (43). This task is arduous. We shall 
employ (31) to circumvent this difficulty. This formula expresses the mass shift in terms of the spectral 
properties of the supersymmetric partner operator 

i^i+ = -^ + ^ + ^sech2x . (45) 

Fortunately, (45) is a differential operator with a Poschl- Teller type potential well whose spectral proper- 
ties are well known. In particular, the spectrum of Ki^ lacks bound states such that presents only 
one bound state, the zero mode. The spatial reflection symmetry of the problem permits us to obtain 
the phase shifts needed in (31) only from the continuous spectrum eigenfunctions: 

V'fc(x) = CP5u.(tanhx) . (46) 

Here, P^{z) stands for Legendre polynomials and C is an integration constant. 

Choosing C = ^^^'^^Jj^^^^^^'^ ^''^ the asymptotic behavior of the eigenfunctions (46) for x — >• — oo 

is: 

, , . x^oo ^ cosh(7r)r[ifc]r[i - z(l + k)]T[\ + i{l - k)] 
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whereas for x — )• oo it becomes 



r[i-i(i + fc)]r[i + i(i-fc)] 

T[-ik]T[l - ik] ' 



ikx 



Consequently the eigenfunctions (46) describe incident waves coming from the left. The reflection 
and transmission scattering amplitudes read 



Pi 



cosh(7r)r[iA;]r[i - i{l + k)]T[l + i(l - k)] 



TrT[-ik] 



(^) 



Til - i{l + k)]T[l + i{l - k)] 

r[-iA:]r[i - ik] 



(^) 



(Ti+ and p\ 



(^) 



pi+ hold because of the symmetry of The 



(1+) 



Moreover, a\ 

eigenvalues cti+ it = e^*^± of the S-matrix provide the phase shifts d'^^^^k) = d''_^^\k) + 6^_'^^\k) 
such that their derivatives give the spectral density: 



dk 



Re [2^P{-ik) - ^(^ - i(l + A;)) - ^(5 + i{l - k))] 



(47) 



Plugging (47) into formula (34), the one-loop mass quantum correction to this kink (44) can be figured 
out: 



11111 
'4 ' 2 4^ ' 2 2^ 



dk 



d5^^+\k) 



dk 



^2 + 1 _ 2 2+^ 2 



The numerical evaluation of the integral in the previous formula gives: ^ dk 



d5(^+){k) 



dk Y '4 

-0.152774. Therefore, the final response to the kink mass quantum shift in this model is 



obtained: 



-0.237985 



4 Construction of scalar field model families sharing the kink mass 
shift 

In references [36, 37] Kumar built one-parametric families of (H-l) dimensional scalar field theory models 
whose kinks give rise to the same second-order kink fluctuation operator Ki-. Evidently, the mass 
quantum shift of all these kinks will be the same, see (27). In this section we shall go one step farther. 
We envisage the identification of one-parametric families of field theory models whose second-order kink 
fluctuation operators Ki^ are different but share the same supersymmetric partner From (31) it 

is possible to compute the common one-loop mass correction to all these kinks in an easy way. We start 
from a specific operator -fCi+, which we assume to be factorized in the usual form: 

= AiA\ where Ai = + Wiix) . 

dx 

The natural question to ask is: does there exist a one-parametric deformation of a superpotential 
Wi{z)dz supplying the same factorization as the reverse Ki^ factorization? The answer is provided 
by a family of the superpotentials such that: 

^i(^^g) = ^i(^) + g^p"^2s,fa)rf^ ' Hi(x) = -y w,{m. (48) 
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We shall restrict the real parameter C to a range such that Wi{x, C) is non-singular. With this definition 
we have a natural relation Wi{x, ±00) = Wi{x), and the zero modes of the family of operators 

Ki- (C) = A{ {C)Ai{C) = + Wi{x,Cf - W[ {x,C) , where Ai (C) = — + Wi {x,C) , (49) 

are: Mx,C) = e" 

One easily checks that for any Wi{x, C) of the form (48) there is a SUSY partner, Ki-^-, common to 
all the members in the Ki- (C) family: 

Ki+ = Ai{C)A\{C) = A,A\ , VC . 

This can be understood from the results of [42]. Taking into account the relations Ai{C)Ki^(C) = 
Ki+Ai{C) and Ki.{C)A{{C) = a{{C)Ki+, one can immediately say that within such a family, any 
two Schrodinger operators are related by the second-order supersymmetry generated by the second-order 
intertwining generators, Y{C,C')Ki-{C') = Ki-{C)Y{C,C'), where 

YiC,C')=Al{C)Ai{C'). (50) 

The idea now is to reconstruct a scalar field model family from this C-parametric family of Schrodinger 
operators. The zero modes are known and consequently the kink solutions of each member of the family 
are obtained through spatial integration: (f)K{x,C) = ipoi^jC) d^. Finally, taking into account that 
4> = 4>k{x, C) are injective functions, they can be inverted, x = in order to derive the potential 

U{(j)) for each family member in the range of the kink profiles: U{4>) = \'4^q[(I>}^ {(!>)]■ We remark that the 
particle masses in all the family members are equal because their kinks fluctuate according to operators 
with the same SUSY partners. 



4.1 A family built from the sine-Gordon model 

The superpotential allowing the factorization of the second-order sG-soliton fluctuation operator (33) is 
Hi (x) = — log cosh X. Thus, the family of superpotentials, all of which share the free particle Hamiltonian 
as the SUSY partner is: 

Si (x, C) = — logfCcoshx + sinhx) , Wi{x,C) = — ; htanhx , (51) 

^ ^ ^ ' ^\ ^ / C + tanhx 

where C is assumed to lead to non-singular VFi(x, C). The second-order soliton fluctuation operator (18) 
arising in each member of the family and their common supersymmetric partner (20) are respectively 

Thus, we have = 1, Vi-{x,C) = (sinh3:+Ccoshx)2 ' = and for later use we list: (Ui_(x,C)) = 

—4, 5^^'^\k) = 0. From (17) and (25) we immediately write the zero mode and the C-dependent kink: 

1 , 2 C tanh % + l 

iI)q{x,C) = - . , (t)K{x,C) = arctan , 

G cosh X + smh x y/C^ — 1 y/C^ — \ 

where the restriction > 1 must be assumed in order to deal with real kinks. In fact, if < 1 
Vi-{x,C) would be singular. Finally, from (26) we also identify the family of potentials 

^ 2C2(C2 _ 1)2 [(^' - 1) (v'C" - 1'^) - 1 sin (Vg^ - 1,/.)] ' . (53) 
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The absolute minima of this function 



^v' = _ ^ {arctan [^/C^-l] + vrj} , j e 



are equivalent in the sense that all of them have the same curvature: v\ = | (j) = 1. As pre- 



"V 



dieted, the particle masses in the whole family are the same. The classical kink energy, however, differs, 
E[(I)k{x,C)] = ) but the computation of the one-loop kink mass shifts by means of formula (31) 
confirms that 



t.E\4>K(C)\ 1 , ,1 



1 r°° 



^(-4) + 2 



1 

TT 



is C-independent and equal to the value found in the sine-Gordon model. Why is this so is revealed by 
consecutively performing a translation and a dilatation in field space. First, the translation 



i = 4> -\ — arctan \fc^ 

\ C — 1 



produces the following modification of the potential: 



U{i, C) = ^(^^ sin2[VC2-l i] . 



Accordingly, the field dilatation ^> = 2\/C^ — 1^ shows that the family of potentials becomes 

U[^,C) = ^^^l_^^ {l-co^^) . (54) 

All of them incorporate the sine-Gordon potential times a C-dependent scale factor. The effect of this 
scale is to deliver the correct C-dependent kink energy by multiplication with the sine-Gordon kink 
energy which is 8. This statement is tantamount to saying that the scale fixes the norm of the zero mode 
because the kink energy density is precisely the square of the zero-mode wave function. 

The requirement > 1 suggests the following reparametrization: cothr = C. In terms of the new 
parameter the superpotential (51) and the function Wi{x^C) become 

Hi(3;, r) = — log — - — cosh(3; -|- r) , WAx^t) = tanhix + t) , 
smhr 

showing that the parameter merely gives rise to a translation in the spatial coordinate of the proper 
sine-Gordon counterparts (modulo a r-dependent integration constant in the case of the superpotential 
H(x, r)): H(x, 0), Wi{x, 0). The kink fluctuation operators Ki^{t) = — ^ -|- 1 — ^^^^^2^^_^^^ are manifestly 
isospectral, the eigenvalues of any i^i_(r) operator are clearly independent of the r parameter. The 
reconstruction of the field theory model from the zero mode, ipo{x^ r) = sech {x+t), and the corresponding 
kink, (j)K{x,T) = 2 arctan tanh , ends in the r-independent potential: U{(I),t) = | cos^ (/> = j(l + 
cos 20). Thus, all the members in the family are completely equivalent to the original sine-Gordon 
model, although a translation in | and a rescaling by 2 are needed. The alternative (and equivalent) 
choice of Hi(x, r) = — log ^jg^ cosh(x -|- r) as superpotential leads to a family of models identical to the 
family proposed by Kumar, discussed above. In summ, translations and dilatations in field space can be 
reinterpreted as translations and dilatations in coordinate space, a situation encompassed in the soldering 
of internal and external variables by topological defects. 

To end this subsection, we mention that in correspondence with (52), the family of isospectral kink 
fluctuation operators (r) can be traced back to the free particle by means of the displaced Darboux 
operators A(t) = ^-j-tanh(x-l-r) and A^(t), as described in references [42, 43, 44]. As a consequence, any 
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pair Ki- (r) and Ki- (r') can be related by the second-order intertwining operator Y{t, t') = A^t)A{t'), 
Y{t,t')Ki-{t') = Ki-(t)Y{t,t') coherently with the discussion of a generic one-parametric family (49). 
The peculiarity of the present family of isospectral kink fluctuation operators is, however, that any two 
of them with r 7^ r' can also be intertwined by the first order operator X{t,t') = ^ — tanh(a; + t) + 
tanh(x-|-r')-|-coth(r— r'), X{t,t')Ki^{t') = Ki^{t)X{t,t'). The product of both, first and second order, 
independent intertwining operators thus generates the non-trivial Lax integral Z{t) = A{t)-^A^ (t) for 
the reflectionless Schrodinger operator Ki_(t): X{t,t')Y{t' ,t) = Z(r) -|-coth(r — r')i^i_, [Z(r),-fCi_] = 
0, see [42, 43] for the details and a discussion of the supersymmetric structure associated with a pair 

i^l„(T), Ki-(t'). 

4.2 The model kinship 

Finally, we shall consider a family of models implicitly constructed by Kumar in [36, 37]. This family is 
based on the (j)^ model whose second-order kink fluctuation operator factorizes in terms of the function 
Wi(x) = 2tanhx. Therefore, application of formula (48) offers the following C-parametric family of 
superpotentials and their derivatives: 

H2(x,C)=log n ; — : : r , W2{x,C) = 2tanh.x 



12Ccosh^x — 2(3sinh3; -|- sinh3x) ' ' 3C — (2 -|- sech^x) tanhx 

The second-order kink fluctuation operators and their supersymmetric partners are 

(P , 6sech^x[4-F9C2 -8sech2x-hsech^x-F6C(-2-Fsech2x)tanhx] 



dx2 [_3C + (2 + sech^x) tanh x] 2 

£■ 

K2+ = J -h4-2sech2x . (55) 

ax^ 

From this we read ^1 = 4 and 

y (J) - 6 sech^x [4 + 9C^ - 8 sech^x + sech^'x + 6C{-2 + sech^x) tanh x] 
2-{x, )- [-3C + (2 + sech^x) tanhx]2 ' 

a family of very complicated potential wells indeed. In Figure 6 the potential well of V2-(x,C) has 
been plotted for several values of C. The condition |C| > | must be imposed to avoid singularities in 
the expression of V2-(x). In this range, we find that (V2_(x)) = —12. The K2+ operator, however, 
is a transparent Poschl- Teller Schrodinger operator that is identical (C-independent) for all the family 
members. The only bound state eigenvalue is = 3, whereas the phase shifts in the eigenfunctions of 
the continuous spectrum are: 6^'^~^\k) = 2arctanp The K2-{C) operators in turn exhibit two bound 
states: the zero mode, (wg = 0) and the SUSY paired state to the iC2+ bound state with energy (the 
same for all the family members) ujf = 3. From formula (17) we identify the zero modes and, implicitly, 
the kink solutions: 

cosli oc f cosh oc 

ipo{x,C) = o , cI)k{x,C) = / dx o . 

12c cosh X — 6sinhx — 2sinh3x J 12C cosh x — 6sinhx — 2sinh3x 

(56) 

Although (56) tells us that the kink solutions exist for any value of C, \C\ > |, the integral in (56) 
cannot in general be expressed in terms of known analytical functions and/or its computation involves 
expressions too complicated to be of practical use. In Figure 6, however, we have depicted the kink 
solutions for several values of the parameter C obtained by means of numerical integration. The value 
C = 6 is exceptional in the sense that the integral (56) is feasible in terms of elementary functions and 
in this case we find the very explicit analytical expression for the kink solution: 



4>k(x, 6) = 

' ^ 960 



/ — — 3 -|- 2 tanhx 

6\/ 15 arctan -— h 10 log[3 cosh x -|- sinh x] — 5 log[5 -|- 7 cosh 2x — 3 sinh 2x] 

V 15 
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Figure 6: Potential wells arising in the second-order kink fluctuation operator K2- (left) and kink solutions (right) 
for several values of C. 



The associated potential energy densities characterizing these field theory models cannot be written 
explicitly because we do not know the analytical form of the kinks. Therefore, we are blind about how 
to invert the kink solution. We do have, however, all the ingredients to compute the one-loop kink mass 
shift in these models. Formula (31) reveals that 



/^E{(l)_ 



1 

4^* 



-12) + 



1 

2^ 



dk 



2^/FT4 
A;2 + l 



-12) + 2 -2 



1 



2^/3 vr 



Thus, all the kinks in this family undergo the same mass shift as the kink of the (/)^ model. 

To finish this sub-Section it is useful to write the above-discussed family of potentials in a way 
that trades the C parameter by spatial displacements of potential wells and kinks. The redefinition 



C 



(tanhr + cothr) is similar but slightly more complicated than the analogous transformation 



performed in the sub-Section 4.1. In the new form the superpotential and the function W2 become 



The zero modes look simpler. 



log 



coshj; 



2tanhx + 



3cosh(x + 2r) + cosh(33; + 2r)_ 
Ssech^x 



cothr + (2 + sech^rE)tanhx + tanhr 



(57) 



coshx 



3cosh(2; + 2r) + cosh(3x + 2r) ' 



dx "00 (^) ''") 



12 ' 



but the potential wells in the family of Hamiltonians 



K2-{t) 



dx^ 



+ 4 



12(3 + 4cosh(2x) + cosh[4(x + r)]) 
(3cosh(a; + 2r) + cosh(33; + 2r))^ 



(58) 



show two minima with a distance between them modulated by the translational parameter r. In fact, 
this Hamiltonian is the ki = 1, K2 = 2, xi = 0, X2 = — r member of the family of two-gap Schrodinger 
operators depending on two scaling and two translation parameters ,to be discussed in the next sub- 
Section. 

In a generic case, in accordance with (50), any pair K2-{t) and K2-{t') of the Schr5dinger 
type operators of the form (58) is intertwined by the second order operator Y{t,t') = A2{t)A2{t'), 
Y{t,t')K2^{t') = K2-{t)Y{t,t'), where A2{t) = ^ + W2{x,t) with W2(x,r) given by equation (57). 
Bearing in mind that the "virtual" system (55) is such that the intertwining relation can be related to the 
free particle Schrodinger operator, the pair of K2-{t) and K2-{t') can also be intertwined by the fifth or- 
der operator X(r,r') = 4(^)4ii^i^2(r'), X{t,t')K2-{t') = K2-{t)X{t,t'), where Ai = ^+tanhx. 
Putting r' = r in the last intertwining relation we see that X{t,t) is nothing else than the nontrivial 
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Lax integral for K2-{t), [X{t,t), K2-{t)] = 0. We refer to [12], where the pecuhar supersymmetric 
structure associated with the isospectral pair of reflectionless Schrodinger operators K2-{t) and K2-{t') 
is discussed in detail. In particular it is shown there that the fifth-order operator X{t, t') can be reduced 
to an independent intertwining operator of differential order 3. 



4.3 From two-gap potentials to scalar field models 

Within the context of supersymmetry, in reference [12], an analysis has been made of the most general 
Hamiltonian, depending on two scaling and two translation parameters, such that the reflection scattering 
amplitude is zero and there are two bound states in its spectrum: 

I 2 2 2^ At|csch^[K2(a:-X2)] + sech^[Ki(x-3;i)] 

J^2- = ^2-^2 = — -r^ + 1^2 — 2{K2 — Ki) -. —. -. ■ — f ttttj , (59) 

ax^ (H2 coth[K2(2; — X2)\ — Ki tanh[K2(2; — a;i)J 

where with no loss of generality we choose K2 > ^- Henceforth, we have ^2 = ^2) (^2-) = — 4(ki + K2), 
and the four-parametric family of potential wells is: 

V2 (x Ki K2 xi X2) = -2iK^ - r?) '^2Csch^[K2(a;-X2)] + Kf sech^[Ki(x-xi)] 
' ' ' ' 21 (/^2 coth[K2(x — a;2)] — Ki tanh[K2(a; — xi)])2 

This second-order operator is factorized by means of the first-order operator A2 = + W2{x) and its 
adjoint, where the superpotential and its derivative are respectively 



l2{x;Kl,K2,Xl,X2) = log 



2 csch[K;2(x — X2) 



K2 coth[K2(x — X2)] — Ki tanh[«;i(x — xi)] 



.... , Kn - Kif^i sech^[Ki(2; - xi)] + K2 coth[K2(x - 2:2)] tanh[Ki(x - xi)]) 

VF2(x;Ki,K2,Xi,X2) = —7— , . r . -(ei) 

K2 coth[K2(x — X2j] — ^1 tanh[Ki(x — xijj 
The ground state, read from the exponential of the superpotential, is a zero mode: 

, / N 2csch[K2(x-X2)] 2 n 

V'o(x;Ki,K2,Xi,X2) = —7 7 7-7 z rr , Wq = • (62) 

K2 coth[K2(x — X2j] — Ki tanh[Ki(x — Xlj] 

The supersymmetric partner operator is surprisingly simple and is related to the family of one-gap 
operators: 

K2+ = A2AI = -— + kI- 2K?sech2[Ki(x - xi)] . (63) 

It is in fact the first member of the hierarchy of reflectionless Poschl- Teller potentials with threshold of 
the continuous spectrum at k^- This operator presents only a bound state '(/'('^2_^2)(x) = sech[Ki(x — xi)], 

with energy — ti\. One easily shows that the first excited state of K2- is paired to this state by 
supersymmetry: 

V'i(x; Ki, K2, xi, X2) = ^2^ech[Ki(x — xi)] , ujI = k^ — k\ 

The total phase shift is also very well known, 5'^'^~^\k) = 2arctan^ and hence 6^'^~\k) is easily obtained 
using supersymmetry. We remark that, unlike the case of one-gap potentials, here there is no a single 
superpartner operator K2+(ki, K2, xi) for all the family operators A'2_(ki, K2, xi, X2) (although in prac- 
tice, the dependence in K2 of K2+ merely shifts the energy levels globally ). Instead, what we find is 
a reduction from a four-parametric family of two-gap operators to a two-parametric family of one-gap 
operators. 



To cope with the opposite inequahty, k-^ > K2, one merely exchanges the roles of ki and /t2. 



26 



Before discussing the scalar field models related to the family of K2-{ki, K2,xi,X2) second- 
order operators it is convenient to analyze some aspects of the family of associated potential wells 
V2-{x, Ki, K2,xi,X2)- In order to elucidate the asymptotic behavior at large values of xi and X2 we 
write the potentials as the sum of two terms: 

V2-{x;Kl,K2,Xl,X2) = V2^\x;Ki,K2,Xi,X2) + V2^{x;Ki,K2,Xi,X2) 

where 

V^_'[x;Ki,K2,Xi,X2) = - 



(coshacosh[Ki(a: — xi)] coth[At2(a; — X2)] — sinhasinh[Ki(x — xi)])^ 



T^(2)/ ^ 

V^_'{x;Ki,K2,Xi,X2) = 



(cosha cosh[K2(x — ^2)] — sinhatanh[Ki(x — xi)] sinh[K2(x — X2))])^ 

The new parameter a is defined in term of the old ones such that: tanha = ^. When K2{x — X2) 3> 
we have that 

V2^\x; Ki, K2, ^1, X2) = — 2K^sech^ ki{x — xi) — sign(x — a;2)Q; 
Analogously, for ki{x — xi) 2> the behavior is: 

V2^{x; Ki, K2, xi,X2) = — 2K|sech^ K2{x — X2) — sign(x — 2:1)0; 

In summ, in this range of parameters the potential wells of the family resemble two onc-boimd state type 
Poschl- Teller potential wells, as if they were describing two separated kinks appearing in two sine-Gordon 
models with different masses: m\ = and = /C2- 

The scalar field theory family is reconstructed from the kink profile identified through spatial integra- 
tion of the zero modes above. It is clear that, like in the = 1 models of previous sections, there is an 
insufficient arsenal of analytical functions to achieve this task without relying on numerical integrations. 
Nevertheless, there are special points in the parameter space {ki, K2,xi,X2) where analytical computation 
of these integrals is possible, leading to highly distinguished scalar field models. 

1. If Ki = 2^2 = 1 and xi = X2 = 0, the kink profile is simply: (t)K{x) = tanhx and we encounter the 
(/)4-modeh U{(t>) = \{<t)'^-\f. 

2. If Ki = ^K2 = c, where c G M is a real constant, and xi = X2 = 0, the kink profile is a dilatation of 
the usual kink in a re-scaled space: ^k{x) = ^ tanh(ca;). Now, the scalar field model responds to 
the potential energy density: U{(f)) = ^{c^4>^ — 1)^, a re-scaled model. 

3. If Ki = \k,2 = 1 and xi = 0, X2 = arctanh g^^ , the kink profile is 

A f ^ R . V5(l + 2e2^) , Vs. 2 + 5e2- + 5e 
<\>k{x) = V3 arctan + — log ^ ^^x^ 



Ax 



4. If Ki = 1, K2 = 4 and xi = X2 = 0, the kink profile is 

(t>K{x) 



1 r. „ /3 . -2 + 3e2^ 



5 -|- 6\/5 arctan — 5 tanh x 

50 L ^ 

This function and the previous one are too complicated to invert. 
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5. At the xi —7- ±00 hmit the kink profile becomes: 

(Pk{x) = , „ „ arctan , • 

The corresponding field model obeys to a potential energy density of the form: 
= 2^ 2 sin^ ^"^ „^ — ^ (/), a sine-Gordon type of scalar field model. 

In order to compute the one-loop mass shifts by using the SUSY generalized DHN formula we now 
pass to the next level of the SUSY hierarchy. To start at a new level we define the operator: 

Ki- = 7^2+ + Kf- kI = — — r + nf- 2k? sech^ [ki(x - xi)] = A\Ai . 

Ki- is again the first transparent Poschl- Teller Hamiltonian with the origin of energies displaced in 
Ki — K2 with respect to K2+. The first-order operators allowing the factorization Ki^ = A[Ai are: 

A\ = — ^ + Wi(x) , Ai = ^ + WUx) where Wi(x) = ki tanhkifx - xi)] (64) 
ax ax 

Clearly vf = and there exists a bound state associated with the eigenvalue of -ftr2+5 which 

is a zero mode of Ki^. The SUSY partner operator in this level is 

Ki+ = AiA\ = -—^+kI 

such that we end in the free particle operator with zero phase shift: J*^^"*") = 0. From (34), the kink mass 
shift for any kink, analytically known or not, in the family can be computed to find: 



1/. , N , 1 n>. ~'2. ------ '^1 



-2{ki + K2) + 2k2 2Ki^k'^ + K^ 



-|- K2) H \ KT) — k\ arcsin ■ 

vr vr V ^ 



As a test, we stress that in the k\ = 1, K2 = 2, xi = 0, X2 = 0, case, which is the (j) -model, the one-loop 
kink mass shift ^^^^'^^ = — ^ is reproduced. We finally mention that this result agrees perfectly 
with the answer obtained using the Cahill-Comtet-Glauber formula: 



^Td ^^2 2 K2 

where the CCG angles are given in terms of the bound state eigenvalues and the continuous spectrum 
threshold. 

To conclude this sub-Section, we note that any isospectral pair K2~ (ki, K2) 2:1, X2) and K2-{ki , K2 ; X]^ , X2) 
of reflectionless Schrodinger operators of the form (59) can be related in a generic case by two indepen- 
dent intertwining operators of differential orders 3 and 2. We do not display their explicit form here 
but simply refer to [ ], where the corresponding associated supersymmetric structure can also be found. 
In a special case, when the relative shift parameters xi — x'l and X2 — x'2 are related by the equality 
Ki coth(xi — x'^) = K2 coth(x2 — X2), the independent irreducible intertwining operators are, instead, the 
operators of orders 4 and 1, 

>4(a;i,X2,x'i,X2) = ^2 t (xi,X2)^l(xi)Ai(Xi)A2(x'i,X2) , (65) 

Xi(xi, X2, x'l , X2) = ^ — h — VF' — Ki coth(xi — x']^) , (66) 
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where W = W2{x] ki, K2, xi, X2) — Wi{x] ki, xi) with W2 and Wi defined, respectively, in (61) and (64), 
while W is defined in the same way as W but with xi, X2 replaced by x'^ and X2. In particular, this 
means that the second order small kink fluctuation operator of the 0^ model with U{(j)) = ^{(j)'^ — 1)^, to 
which the choice ki = = 1, xi = X2 = does correspond, can be related by exotic supersymmetry 
involving the supercharges of orders 4 and 1, composed from (65) and (66), with the kink fluctuation 
operator of the form (59) with ki = = 1 and shift parameter xi given in terms of the arbitrary shift 
parameter X2 by a relation xi = arccotanh(2cothx2). 

5 Conclusions and further comments 

The main conclusion to be drawn here is the fact that the intrinsic supersymmetry attached to the 
stability operator of BPS kinks is helpful in the task of evaluating the one-loop mass shifts induced by 
kink fluctuations. Throughout this work we have applied this idea in increasing order of difficulty: first, 
we addressed very well known models from this new strategic point of view to test the method and gain 
confidence in its effectiveness. Second, we explored ignotum models where the new procedure showed 
the answer, impossible to reach by means of traditional approaches. Finally, within this new framework 
we have analyzed families of models with a very rich supersymmetric structure. Again, the results were 
extremely satisfactory. 

In hindsight, these developments lend us to speculate that the extension of the method to scalar field 
models with more than one scalar field will pay even better dividends. Even though the heat kernel 
expansion allowed the computation of BPS kink mass shifts in some interesting models, [24, 31], there 
are other more interesting two-component BPS Wess-Zumino kinks, [45, 46, 47] where this computation 
has not yet been achieved within the purely bosonic framework. Starting with the models where some 
results are available we plan to address the computation of Wess-Zumino kink mass shifts using the SUSY 
quantum mechanical methods developed in this paper. 

Finally, we look forward to treating BPS topological defects in higher dimensions similarly, e.g. 
the BPS self-dual vortices in the = 2-dimensional Abelian Higgs model at critical coupling between 
Type I and II superconductivity. To accomplish this task, one must deal with supersymmetric quantum 
mechanical systems in A^ dimensions. The spectral problems are thus much more involved but a fairly 
good quantity of information is available, see e.g. [ JS] to see a recent review. 
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